Magnetic Susceptibility of the Kagome Antiferromagnet ZnCu 3 (OH) 6 Cl2 
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We analyze the experimental data for the magnetic susceptibility of the material ZnCu3(OH)eCl2 
in terms of the Kagome Lattice Heisenberg model (KLHM), discussing possible role of impurity 
spins, dilution, exchange anisotropy, and both out-of-plane and in-plane Dzyaloshinsky-Moriya 
(DM) anisotropies, with explicit theoretical calculations using the Numerical Linked Cluster (NLC) 
method and exact diagonalization (ED). The high-temperature experimental data are well described 
by the pure Heisenberg model with J = 170-K". We show that the sudden upturn in the susceptibility 
around T = 75K is due to DM interactions. We also observe that at intermediate temperatures, 
below T = J, our calculated susceptibility for KLHM fits well with a power law J 1-0 - 25 . 

PACS numbers: 75.10. Jm,05.50.+q,05.70.-a 
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Frustrated magnetic systems represent outstanding 
challenges in condensed matter physics While much 
progress has been made theoretically over the last few 
decades 0], key issues related to concrete lattice mod- 
els remain unresolved. Computational methods are re- 
stricted to small system sizes or high temperatures, 
whereas the starting point for field theory methods, 
which maybe asymptotically exact at low temperatures, 
involve drastic approximations, which makes their valid- 
ity for realistic lattice models unclear. Experiments can 
bridge this gap by providing information all the way from 
high to low temperatures. 

The spin-half Kagome-Lattice Heisenberg Model 
(KLHM), is one of the best studied frustrated quantum 
spin models. Much of our knowledge of the system comes 
from studies of finite-size periodic clusters. Exact Diago- 
nalization studies 0] suggest that the model has a short 
spin-spin correlation length, and a spin gap of approx- 
imately J/20. In the finite systems, there are a large 
number of singlet states below the lowest triplet and their 
number grows with the system size Q ■ Various Resonat- 
ing Valence Bond (RVB) scenarios have been invoked, 
and the finite size spectra have also been interpreted to 
imply deconfined spin- half excitations. 

In light of these theoretical results, the experimental 
behavior of the recently synthesized Kagome-Lattice ma- 
terials ZnCu3(OH) 6 Cl2 are highly unexpected @,@|. The 
high temperature inverse susceptibility data was found 
to obey a Curie- Weiss law, with a Curie- Weiss constant 
of about 3Q0K. Yet, there appears no spin-gap either 
in the susceptibility or the specific heat or the neutron 
spectra or in the nuclear spin-lattice relaxation T\ , down 
to temperatures below lOOmiV At the lowest temper- 
atures the susceptibility saturates and the specific heat 
shows power-law behavior in temperature. The latter is 
suppressed by magnetic fields, showing it to be magnetic 
in origin. A possible interpretation of the substantial 
rise in susceptibility at low temperatures is that it is due 
to impurity spins outside the Kagome planes, caused by 



substitutions of non-magnetic Zn sites with Cu [7]. How- 
ever, the fact that the muon shift K tracks the bulk sus- 
ceptibility Xi implies that the susceptibility is not due to 
inhomogeneities localized away from the Kagome planes, 
but rather a bulk behavior of the system. 

Here, we use the triangle-based NLC method [1] and 
ED, to calculate the magnetic susceptibility of the KLHM 
including various perturbations such as dilution, ex- 
change anisotropy and Dzyaloshinsky-Moriya anisotropy. 
We find that the high temperature susceptibility is indeed 
in excellent agreement with that of KLHM. However, the 
Curie- Weiss fit is fortuitous. The true Curie- Weiss be- 
havior (with Tew = J) is only valid for T > 10J, a re- 
gion inaccessible to experiments. When Curie- Weiss fits 
are made to data for T < 2 J, the effective Curie- Weiss 
constant increases in units of J, and hence the true J is 
smaller than reported in Ref. [a, Q . The high tempera- 
ture data fits well with J = 170K. 

The sharp upturn in the susceptibility around 75K 
cannot be accounted for by the KLHM. We find that 
only a Dzyaloshinsky-Moriya (DM) anisotropy can ac- 
count for the sharp upturn in the experimental suscepti- 
bility. A way to distinguish this from an impurity spin 
would be by the anisotropy. The excess susceptibility due 
to DM would be strongly anisotropic, and this anisotropy 
is temperature dependent. 

The NLC method allows us to calculate the suscep- 
tibility for the KLHM accurately down to T = 0.3J, 
without extrapolations. In contrast, the high temper- 
ature expansions (HTE) fail to converge below T = J 
. Some extrapolations of NLC and HTE suggest that 
shortlybelow T = 0.3 J, the susceptibility begins to turn 
down 8, 9]. Unfortunately, these extrapolations to lower 
temperatures become subjective, and hence we focus only 
on the region T > 0.3 J. All NLC calculations reported 
below are done using the triangle based expansion, where 
contributions to the thermodynamic system from graphs 
with complete triangles are included [8j. As appropriate, 
two highest order calculations are shown. 



2 



The Heisenberg Hamiltonian in field h is 

H = jY^a i a j -gn B h y £ i sj. 

(«.i> i 



(1) 



In our calculations we set J = 1, Boltzmann constant 
k = 1 and g/i^ = 1. The susceptibility per spin is 



T <9 2 lnZ 



X N 



dh 2 



(2) 



h=0 



where Z is the partition function. The molar suscepti- 
bility, measured experimentally, is related to our suscep- 
tibility per spin by the relation, \moiar = C\, where the 
constant C = NAg 2 fJ? B /kJ — 0.3752g 2 / J in cgs units. At 
high temperatures the susceptibility has a Curie- Weiss 
form with Curie- Weiss temperature T cw — J/k. As dis- 
cussed by Zheng et al. [l(J, asymptotic Curie- Weiss be- 
havior is only valid for T ^> J. At temperatures of or- 
der J, one can define an effective temperature dependent 
Curie- Weiss constant as 



t; 



eff 



-T 
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dx/dT 
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If one was to fit to a linear curve in the vicinity of 
some temperature T and use the intercept to estimate the 
Curie- Weiss constant, one would get Tj~*. This quantity 
for KLHM is shown in Fig. 1. Note that for T < 2 J one 
has T c c 5 w 2J. This means that J for ZnCu 3 (OH) 6 Cl 2 
is smaller than estimated in Refs. 0, @| . In the inset of 
Fig. 1, we show a comparison of the experimental sus- 
ceptibility with the theoretical one with J = 170K. The 
agreement is excellent at high temperatures. 




FIG. 1: Temperature dependent effective Curie-Weiss tem- 
perature Tel, for the KLHM. The inset shows fits to the 
experimental data of Ref. [f|. For both experimental cases 
we have taken J = 170A'. The g-factors needed for the fit 
are 2.19 for the data from Ref. [5j and 2.33 for the data from 
Ref. [6|| . On the theoretical side, the two highest order NLC 
results are shown. 

While the theoretical susceptibility for the KLHM has 
a weak upturn at T near 0.3J, it is inconsistent with 



the pronounced rise in the measured susceptibility, which 
grows significantly below 75K and only saturates to a 
much higher value at temperatures below IK. The ex- 
perimental behavior could arise from several sources. 
First, there maybe external impurity spins, such as those 
caused by substitution of Zn sites with Cu pj], which 
are only very weakly coupled to the electronic spins in 
the Kagome planes. If present, they would add a Curie 
like 1/T term to the susceptibility at high temperatures. 
Fig. 2 shows that considering 3 percent impurities allows 
us to extend the fit down to T = 0.3J, but there re- 
mains a sharp rise at lower temperatures. A 6 percent 
impurity concentration eliminates the upturn but makes 
the data inconsistent with the KLHM. Hence, free im- 
purity spins by themselves, cannot account for the rise 
in the data. Besides, as argued by Ofer et al. [Bj], in 
this scenario, it would be difficult to understand why the 
muon-shift tracks the bulk susceptibility. One would ex- 
pect local probes to show inhomogeneous behavior and 
deviate from the bulk averages. 
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FIG. 2: Comparison of KLHM susceptibility with experimen- 
tal susceptibility after subtracting Curie contributions from 3 
and 6 percent impurities. 

A second potential reason for the sudden ferromag- 
netic behavior could be dilution due to substitution of 
Cu sites in the Kagome-planes by Zn. The missing spins 
on the lattice could create local moments in the singlet 
background and cause a Curie-like susceptibility to arise 
at low temperatures. They could weakly couple to each 
other and hence saturate at some much lower temper- 
ature. In order to investigate this possibility, we have 
studied the KLHM with quenched dilution. We assume 
that at each site, we have a hole with probability c and 
a spin with probability 1 — c. The holes are fixed in their 
position and extensive quantities are averaged over the 
random configurations C, using the relation 



(0)=Y / P(C)0(C) 



(4) 



where P(C) is the probability of the configuration C. 
The susceptibility for several values of hole concentra- 
tion is shown in Fig. 3. We note that at these interme- 
diate temperatures, holes simply lower the susceptibility 
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and do not lead to any ferromagnetic tendencies. Thus 
the static holes do not provide an explanation for the 
experimental data, either. We have also calculated the 
susceptibility of KLHM with exchange anisotropy of Ising 
and XY type. Neither can account for the experimental 
behavior. 
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FIG. 3: NLC results for the magnetic susceptibility of the 
KLHM with hole concentration c ranging from to 10 percent. 

We now turn to the Dzyaloshinsky-Moriya (DM) inter- 
action. In Kagome magnets (which we assume lies in the 
x — y plane), both out of plan e (D z ) and in plane (D p ) 
DM terms are allowed fTTL 1T^] . 

xSj) 



H 



DM 



= 5>,(Si 

(id) 



D p .(SiXSj), (5) 



The D z term alternates between the up and down point- 
ing triangles of the Kagome lattice. Its sign can be set 
by demanding that for the up pointing triangle shown 
in Fig. 4b, a positive D z multiplies (Si x S 2 ) 2 - The in- 
plane DM term D p is perpendicular to the bonds and 
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points inward towards the center of the triangles 
Since the DM terms break spin rotational symmetry, we 
need to calculate separately the susceptibility with field 
along z (xz) and in the x — y plane (x P ). The powder 
susceptibility Xa is given by Xa = \{^Xp + Xz)- 

With the DM terms, the different S z sectors are cou- 
pled so we are able to do NLC calculations only up to 
6 triangles [8]. Unfortunately, the convergence is poor, 
and hence we turn to ED of clusters with 12 and 15 sites 
(with periodic boundary conditions) [§] to study the ef- 
fects of DM interactions. We find that a pure D z term 
supresses the in-plane and z susceptibilities with respect 
to KLHM. On the other hand, a pure D p enhances both 
in-plane and z susceptibilities with respect to KLHM. 
Hence, the competition between D z and D p can produce 
varying results. Once D p ^ the sign of D z is also rele- 
vant 
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In all cases studied, the DM terms produce z 
susceptibilities which are larger than the in-plane ones. 

As long as D p > \D Z \, we find an enhancement of the 
powder susceptibility with respect to the KLHM. As an 
example, we show in Figs. 4(a) and 4(b) results for the 
powder susceptibilities and anisotropy when D p = 0.3J 



and D z — — 0.15J, — 0.3J, respectively. Note that in 
the theoretical calculations the g-factor is assumed to be 
isotropic. Since the z-susceptibility rises very rapidly, an 
anisotropic ^-factor enhanced along z will cause an even 
more rapid rise and will lead to agreement with experi- 
ments with a smaller DM anisotropy 

Our main point with Fig. 4 is to show that even though 
more information from the material ZnCu3(OH)gCl2 may 
be required in order to make a precise estimate of the DM 
anisotropy, DM interactions are essential to understand- 
ing that material. We predict that single crystal mea- 
surements should see an upturn in the anisotropy when 
the powder susceptibility departs from the KLHM re- 
sult. Such a behavior has already been seen for spin-5/2 
Kagome system KFe 3 (OH) 6 (S0 4 )2 0. 
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FIG. 4: (a) 15 sites ED calculation of the powder suscepti- 
bility with a DM term D p — 0.3 J and different values of D z , 
compared with experiments and pure KLHM. (b) Anisotropy 
vs temperature for D p = 0.3 J and different values of D z . 

Here, we cannot address the measurements at T <C J. 
An interesting question to ask is whether the lack of spin- 
gap seen in the experiments is due to DM anisotropy. In 
magnetically ordered systems, the DM term often creates 
a spin-gap in the spin-wave spectra. The KLHM is an 
interesting system where the spin-gap maybe non-zero 
but singlet excitations maybe gapless 0,0). In that case, 
the DM term can cause the spin-gap to vanish by mixing 
the singlet states with the spinful states. Alternatively, 
the experimental data can be interpreted to imply that 
the KLHM has no spin-gap in the thermodynamic limit. 

We now turn our attention back to the pure KLHM. 
The log-log plot in Fig. 5 makes clear that around T = J, 
where HTE convergence fails, there is a crossover in 
the temperature dependence of the susceptibility. Be- 
low T = J, the susceptibility is remarkably well fit by 
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a power law in temperature of the form T~ - 25 . This 
non-trivial crossover behavior between the high temper- 
ature Curie- Weiss regime and the very low temperature 
RVB regime deserves further theoretical attention. Ex- 
perimental measurements of wavevector dependent sus- 
ceptibilities, at temperatures below J can shed further 
light on this crossover. 
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FIG. 5: Log-log plot of the susceptibility of the KLHM. In the 
inset we show log-log plots of the entropy of the pure KLHM 
and in the presence of DM terms. Straight lines following 



the data are \ ^ 
5 ~ T 0A (inset). 



(low T), % = V(4T) (high T), and 



A rough power law is also seen in the entropy function 
of the KLHM, which goes as S ~ T aA (Inset in Fig. 5). 
This may help explain the failure of naive Pade extrap- 
olations for the specific heat of the KLHM calculated by 
HTE [9( . It was found that Pade approximations lead to 
specific heat curves, which when integrated give a finite 
ground state entropy. It was suggested that this meant 
the KLHM specific heat has two peaks as a function of 
temperature as discussed by Elser earlier 17j, and HTE 
simply fails to recover the second peak. However, there is 
no strong evidence to support the second peak. Our re- 
sults suggest that instead of a second peak, there maybe a 
non-trivial power-law regime at low or intermediate tem- 
peratures which simple Pade extrapolations of HTE fail 
to capture. As shown in the inset in Fig. 5, the addi- 
tion of DM terms to the KLHM produces an increase 
in the exponent of the entropy from 0.4, bringing it in 
rather good agreement with low temperature measure- 
ments [6|. If the magnetic specific heat for the material 
ZnCu3(OH)eCl2 can be extracted at higher T, the above 
issues can be further addressed. 

In conclusion, in this paper we have studied the mag- 
netic susceptibility for the material ZnCu3(OH) 6 Cl2 by 
comparing it with theoretical calculations for KLHM. 
We find that at high temperatures the susceptibility 
data is well described by KLHM with J = 170K. The 
sudden upturn in the experimental susceptibility below 
T = 75K is not consistent with KLHM but can be ex- 



plained by Dzyaloshinsky-Moriya interactions. Further 
experiments, can verify or refute the existence of such 
terms by measuring the anisotropies in the susceptibil- 
ity. We have also shown that the susceptibility of KLHM 
has an intermediate temperature regime where x goes 
as T~ - 25 . The combined experimental and numerical 
results raise a fundamental question about whether the 
KLHM has a spin-gap, or whether it has a gapless spin- 
liquid ground state [7J, LL8( . 
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Note Added. — Recent NMR experiments by Imai et al. 
[lit find that parts of the NMR spectra show a temper- 
ature dependence like the bulk susceptibility that grows 
rapidly as T is lowered while other parts flatten out and 
then decrease at still lower T. This is exactly what we 
find for Xz and Xv respectively. Thus a simple explana- 
tion for these experiments is that they arise from powders 
oriented perpendicular and parallel to the field. 
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